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ABSTRACT 

Axisymmetric steady-state weakly ionized Hall-MHD Keplerian thin disks are inves- 
tigated by using asymptotic expansions in the small disk aspect ratio e. The model 
incorporates the azimuthal and poloidal components of the magnetic fields in the 
leading order in e. The disk structure is described by an appropriate Grad-Shafranov 
equation for the poloidal flux function -0 that involves two arbitrary functions of for 
the toroidal and poloidal currents. The flux function is symmetric about the midplane 
and satisfles certain boundary conditions at the near-horizontal disk edges. The bound- 
ary conditions model the combined effect of the primordial as well as the dipole-like 
magnetic flelds. An analytical solution for the Hall equilibrium is achieved by further 
expanding the relevant equations in an additional small parameter S that is inversely 
proportional to the Hall parameter. It is thus found that the Hall equilibrium disks 
fall into two types: Keplerian disks with (i) small {Rd ~ (5°) and (ii) large {Rd ^ S"'', 
k > 0) radius of the disk. The numerical examples that are presented demonstrate 
the richness and great variety of magnetic and density conflgurations that may be 
achieved under the Hall-MHD equilibrium. 
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1 INTRODUCTION 

In the present study the Hall equilibrium of thin Keplerian disks embedded in 3D axially-symmetric magnetic fields is 
investigated. Protoplanetary disks have been the subject of numerous investigations in search of possible instabilities that 
may play a significant role in such yet not totally understood phenomena as the outwards transfer of angular momentum, 
and planet formation. However, in most astrophysics-related hydromagnetic stability studies which include accretion, jet, and 
wind effects the models for the underlying equilibrium states are based on the local app roximat i on or a 'cylindrical disk 



that i s uniform in the a xial d irection. A more realistic model of thin disks is adopte d in Reeev ( 19831 ). iKluzniak fc Kita 



I 2000l ). lUmurhan et~all l|2006l ) or equilibrium nonmagnetized rotating disks (see also lOgilviel l|l997h for Keplerian disks of 
ideal MHD plasmas) and is based on an asymptotic approach in the small aspect ratio of the disk. Traditionally, a big variety 
of geometrical configurations of axially symmetric steady-state equilibria are described through a few arbitrary functions of 
the poloidal magnetic flux. The latter satisfies the Grad-Shafr a nov (GS) type equations within bo th one-fluid and multifluid 
m odels of the plasmas (e .g. McClements fc Thvagaraia ( 2001 ): Thvagaraia fc McClements ( 20061 ) . see also the earlier study 
bv lLovelace et al. I l| 19861 ). where in particular thin rotating disks have been discussed and references therein). Such solutions 
provide a base for their stability studies. 

Although the magnetic field configuration in real disks is largely unknown and, as observations and numerical simulations 
indicate, toroidal magnetic component may be of the same order of magnitude, or some times domina t e the p oloidal field due to 
the di ff erential rotati on of the disk (see e.g. iTerguem fc Papaloizou I II 19961 ): iPapaloizou fc Tergueml ([1997|) ; |Hawj«y_&_j<rolik 



( 2002l ): |Proga1(|2003l )). As shown in the present study there are three kinds of possible equilibrium in thin rotating disks which 
can not be reduced from one to the other: (i) pure poloidal, (ii) pure toroidal and (iii) mixed toroidal and poloidal magnetic 
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field equilibria. The first case has been considered within the ideal MHD model by Lovelace et al. ( 19861 ). and Ogilvie (Il997 ) 
The s econd case, i.e. the pure toroidal magnetic field has been investigated within the Hall MHD model bv lShtemler et al 



(2007]). The present study is aimed to describe numerous axially-symmetric equilibria in 3D axially-symmetric magnetic fields 
that contain both toroidal and poloidal components, i.e. case (iii) within the Hall MHD model in thin disk approximation. 

While most of recent research has focused on the magnetohydrodynamic (MHD) description of magneto-rotational insta- 
bility (MRI), the importance of the Hall electric field to such astrop h ysical objects as protoplanetary d i sks has been pointed 
out by Wardle (1999 ) . Since then, such works as 



Balbus 



sucn astrop n ysicai oD.iects as protoplanetary g i sKs nas Deen pomtea 
fc TerauemI j200ll).ISalmeron fc Wardi3 (l2003l).ISalmeron fc Wardlel 



( 2005l ). lDeschl (|2004l ') . [Urpin fc Rudigei] (|2005l ). iRudiger fc Kitchatinovl l|2005l ). and lPandev fc Wardlel (|2008l ) have shed hght 
on the role of the Hall effect in the modification of the MRIs. However, in addition to modifying the MRIs, the Hall electro- 
motive force gives rise to a new family of instabilities in a density-stratified environment. That new family of instabilities is 
characterized by the Hall parameter tth of the order of unity ( tth is the ratio of the Hall drift velocity Vhd to the charac- 
teristic velocity, where Vhd = li^i/H^, U a nd Q,i are the inertial length and Larmor frequency o f the ions, respectiv e ly, an d 
Hi, is the density inhomogeneity length fsee iHubal l|l99ll )V As shown bv lLiverts fc Mondl (|2004 ). and lKolberg et al.l (|2005l ), 
the Hall electric field combined with radial density stratification gives rise to two modes of wave propagation. The first one is 
the stable fast magnetic penetration mode, while the second is a slow quasi-electrostatic mode that may become unstable for 
a short-enough density inhomogeneity length H^- The latter is termed the Hall instability. Density stratifications also play a 
crucial role in the recently discussed starto-rotational instability. That instability is excited by the combined effect of vertical 
and radial density stratification. 

The aim of the present work is, therefore to describe the various Hall equilibrium configurations of magnetic field and 
density in thin magnetized Keplerian disks as a first step towards a further study of their various density stratification- 
dependent instabilities. The following simplifying assumptions are adopted throughout the present work: the plasma is treated 
in a 3D axially symmetric magnetic field within the Hall MHD model. The problem inside the disk is treated within the thin- 
disk approximation, while the effects of the outer region is modeled by the boundary condition f or the p o loidal magnetic 
flux. Influence of the accretion, jet and wind effects are neglected in the thin-disk approximation (jOgilvid ([1997|)). A near 
Keplerian region of the rotating disk is considered that starts at some radial distance away from a central body. The influence 
of the latter on the Keplerian portion of the disk is modeled by a dipole-like contribution (singular at the disk axis) to the 
boundary condition for the poloidal magnetic flux at the near-horizontal disk edges. In addition, axial electrical currents that 
are localized in the central non-Keplerian regions of the disk are taken into account through the resulting toroidal magnetic 
fleld components within the rotating disk. Influence of the interstellar (primordial) magnetic field is modeled by specifying 
the contribution of the latter on the near-horizontal disk edges. 

The paper is organized as follows. The dimensional governing equations, the normalization procedure, and the resulting 
non-dimensional system are presented in the next Section. Section 3 contains an asymptotic analysis of the Hall-equilibrium 
in the small aspect ratio as well as the derivation the appropriate GS equation for the fiux function and the corresponding 
boundary conditions at the disk edge. In section 4 an analytic model is developed for the Hall equilibrium with the aid of 
asymptotic expansions in a newly defined small parameter inversely proportional to the Hall parameter. Numerical examples 
for different regimes of the equilibrium are presented in section 5. Summary and discussion are in section 6. 



2 THE PHYSICAL MODEL FOR THIN KEPLERIAN DISKS 

A rotating thin disks is considered that is under the influence of a general 3D axially-symmetric magnetic fleld, as well as the 
gravitational fleld due to a massive central body. Electron inertia and pressure, displacement current, viscosity and radiation 
effects are neglected in the electrons momentum equation. Consequently the latter is reduced to the generalized Ohm law that 
takes into account the Hall effect. In addition the plasma is assumed to be quasi-neutral. 

2.1 The basic equations 

Under the assumptions mentioned above the physical model of the rotating magnetized disk is given by: 
^-fV.(nV)=0, (1) 
DY 1 

miW— — = -VP-I- -j X B - minV$, (2) 
Dt c 

^^{Pn-'')^0, (3) 

^ + cV X E = 0, V • B = 0, (4) 

at 
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E = — V X B H ^- , j = — V X B. (5) 

c ec Tie 47r 

Here the standard cylindrical coordinates {r, 9, z} are adopted throughout the paper with the associated unit basic vectors 
{ir,ie,iz}; V is the plasma velocity; t is time; D/Dt — d/dt + (V • V) is the material derivative; = —GMc/R is the 

gravitational potential of the central object; — r'^ + z'^; G is the gravitational constant; Mc is the total mass of the central 
object. The electric field E is described by the generalized Ohm's law which is derived from the momentum equation for the 
electrons fluid by neglecting the electrons inertia and pressure; B is the magnetic field, j current density; 7 is the polytrophic 
coefficient, (7 — 5/3 in the adiabatic case). P = Pe + Pi + Pn is the total plasma pressure; c is speed of light; Pi and mi are 
the species pressures and masses {I — e,i,n ), subscripts e, i and n denote the electrons, ions and neutrals, respectively; e is 
the electron charge; rrii = ZrUp, nip is the proton mass {Z = 1 for simplicity). Since the plasma is assumed to be quasi-neutral 
and partially ionized with small ionization degree and strongly coupled ions and neutrals 

He ^ Ui ann,n ^ nn,a — Ue/ [rie + Un) « 1, Vi ~ Vn. (6) 

In real protoplanetary disks, the electron density is determined by ionization versus recombination, and the fractional 
ionization may vary significantly in space. In such cases, rate equations that describe the ionization degree should be employed. 
However, in the present study the ionization fraction is assumed to be constant. This simplification is made to avoid the widely 
uncertain physics of ionization and recombination. Nevertheless, such ap proximation allows on e to roughly estimate the real 



characteristics of the equilibrium system (see a discussion in section 5 in lShtemler et al.l l|2007f )) 



2.2 Scaling procedure 

The physical variables are now transformed into nondimensional variables: 

/nd = (7) 

where / and /„d stand for any of the physical dimensional and non-dimensional variables, while the characteristic scales /* 
are defined as follows: 

v. = fl^r,, t, = $. = V;^ m, = rui, n, = n„, P. = K{m^n^y, j, = — — , E-^ = Yl^, (8) 
Si* Atv r, c 

Here fi* = {GMc/r^Y^"^ is the Keplerian angular velocity of the fluid at the characteristic radius r, that belongs to the 
Keplerian portion of the disk; K is the dimensional constant in the steady-state polytropic law P — Kn' . The characteristic 
values of the electric current and fleld, j, and iJ, , have been chosen consistently with Maxwell equations. The characteristic 
dimensional magnetic fleld _B* is specified below. Note that a preferred direction is tacitly defined here, namely, the positive 
direction of the z axis is chosen according to positive Keplerian rotation. 

The resulting dimensionless system (omitting the subscript 'nd! from the dimensionless variables) is given by: 

^ + V(nV)=0, (9) 

D{Pn- 
Dt 

^ + VxE = 0, V-B = 0, (12) 



Dt Ml ^ (3 Ml' -V',-; (^2 + 22)1/2 

= 0, (11) 



E= -V X B + Hj/i^, j = VxB. (13) 
n 

Here Ms and 13 axe the Mach number and plasma beta, and YIh is the Hall coefficient in the generalized Ohm' law (13): 



cs* ~ \J P, / {rn,n») is the characteristic sound velocity, U — c/cupi and fii = eB, / {m*c) are the inertial length and the Larmor 
frequency of ions, respectively, copi = ^ Aire^an* /nit is the plasma frequency of the ions. The increasing importance of the 
Hall term for weakly ionized disks is apparent as Hh is inversely proportional to the ionization degree a. 

A common property of thin Keplerian disks is their highly compressible motion with large Mach numbers Ms- Further- 
more, the characteristic effective semi-thickness H, = H{r,) of the disk {H = H{r) is the local disk height) is such that the 
disk aspect ratio e equals the inverse Mach number: 
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Table 1. Estimation of the physical conditions at two radii within the disk 



r. AU 


a 


li Km 




ujpi sec ^ 


Qf sec ^ 


1 


1.12 X 10"" 


1.24 X 10^ 


1.03 


241.6 


2.00 X 10"'' 


20 


4.14 X 10"" 


1.25 X lO* 


11.2 


23.82 


2.23 X 10"*' 



T^ = ^-^<1. (15) 

Ms r* 

The smallness of e means that dimensional axial coordinate 1^1 ^ iif is also small, i.e. z/r, ~ e, and the following rescaled 
values of the order of unity in e may be introduced in order to further apply the asymptotic expansions in e (Shtemler et al 
2007): 

C=-, h= , TTH = , (16) 

e e e 

where as was mentioned in the introduction, the rescaled Hall parameter tth is determined in terms of the characteristic disk 
thickness _ff, = er, as the density inhomogeneity length. 



2.3 Physical conditions in protoplanetary disks 



Before turning to the detailed description of the Hall equilibrium of rotating Keplerian disks it is instructive to estimate some 
of the parameters that have been introduced in the previous subsection. Typical protoplanetary disks extend up to the order 
of IQQ AU and consist of molecular gas with characteristic ion masses between SOrrip to 40mp. The temperature within a radius 
of about Q.IAU may exceed lOpK were as further away from the central star the temperature decreases and may reach values 
of about IQK in the outer regions of the disk. As a result, while in the inner regions of the disk thermal collisions provide the 
dominant ionization mechanism, in region s beyond about Q.IAU the onl y sources o f ioni z ation are nonthermal, li ke cosmic 
rays and the decay of radioactive elements l Umebavashi fc Nakano ( 19881 ). Gammie ( 19961 ). Igea fc Glassgold ( 19991 )). Taking 
that into account, the frequently used model is employed, in which the radial temperature profile is given by 



r(r) = 280(: 



-1/2 



■lAU' 

while the column mass density is given by the minimum mass model as 
r .-3/2 



E(r) = 1700( 



lAU' 



g cm 



(17) 



(18) 



The thickness of the disk as well as the neutral number density radial profiles may be calculated from the two functions above. 
The ionization fraction is given by 



(19) 




where 5, the ionization rate due to the non thermal processes, is given by 

C = 10"^^exp(-E/192) + 7 X 10"^^ sec"\ 

and j3r is the dissociative recombination rate given by 

Pr = 1.1 X 10"^(— ^) 



300 A" 



3 -1 

cm sec 



(20) 



(21) 



A det ai led and critical study o f the model as well a s some other suggestions may be found in iHavashi et al. 

I l|l985h . iGammid 

l l99d ). lFromang et all (|2002| ). ISano fc Stonj (|2002l ). Here however, Eqs. (fT7)l - (pT]) are employed in order to obta n a rough 
estimate for the physical conditions in the disk. To do that, attention is focused on two representative radii, r\ — lAU and 
r2 — 20AU. A characteristic magnetic field of O.IG is assumed while the mass of the central star has been assumed to be one 
solar mass. The results are summarized in Table (1). The increase in the rescaled Hall parameter is due to the decrease in 
the neutrals number density that compensates for the increase in the ionization fraction and the angular momentum per unit 
mass. In particular, it is evident that the entire disk resides within the Hall regime. 



3 HALL EQUILIBRIUM IN THIN DISKS IN 3D AXISYMMETRICAL MAGNETIC FIELDS 



In magnetized disks the solution of the steady-state equilibrium problem {d/dt 



portion of a disk by asymptotic expansions in small e with the aid of Eqs. (f9t- (|16|l (similar to lRegev 



0) may be obtained with i n the Kepleria n 
19831 ), lOgil^ (|l997l ). 
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Kluzniak fc Kital l|20od) . IUmurhan et all (|2006l ). and lshtemler et alj (|2007l )). This yields for the gravitational potential in the 
Keplerian portion of the disk: 



(22) 



1 e^C 

To leading order in e within the thin disk approximation, the toroidal velocity Ve is described by the Keplerian law, which 
follows from the leading order radial component of the momentum equation: 

Y^Vgig, Vg = Q.{r)r, n(r-) = r-='''^ Vg 



e , Vr 



(23) 



The asymptotic orders in e in equilibrium relations (|23|l are similar to those in iQgilvid (|l997l ). They may be inferred by 
observing that due to Eq. ((9} (V • (nV) = 0) the axisymmetric fluid momentum per unit volume is expressed in terms of a 
scalar stream function x('", C) 



nV = nVgig + -Vx X \g 



d 



1 d 



)■ 



(24) 



Assuming now that the toroidal velocity component is dominant, and noticing that d/dz ~ imply that x ~ which 
immediately results in the ordering presented in Eq. (|23p (see also discussion just after Eq. (|41|l ). 

By a similar way the divergent free axisymmetric magnetic field B, (V • B = 0) is written in terms of a poloidal flux 
function '^{r,C,) in the following way: 

B = Beie + iv* X is, *(r, C) = *o(r) + e^(r, C), Br = --^, B. = -^+e-^. (25) 
r r oQ r ar r or 

Here the poloidal flux function '3/ is presented as sum of 'i'o{r) and etp{r,Q, which are scaled in e by such a way to produce 
Bz ^ Br ^ to leading order in the thin disk approximation. 

Proceeding further it is noticed that Eq. (|12p means that the electric field is derived from a potential (j){r, C,) 



E = e\/(j}. 



(26) 



Taking now the dot product of both sides of Eq. (|13p with the magnetic field reveals that B • V(/f> = 0, and therefore is a 
function of the magnetic fiux, i.e. 



(27) 



which immediately, due to axisymmetry means that the toroidal component of the electric field vanishes, i.e.. Eg = 0. In order 
to impose that condition. Ampere's law is first employed in order to write the electric current density in the following way: 



. .. 1^, . . IdBg . 1 d /d*o\ 119V d (ldi^\ 

i=Jgig + -SI{rBg)xie, jr ^ tt— , Je = -r[~^l ' 

r € oC r ar \ ar / € r OQ'^ 



dr\rdr'^ r dr 



(28) 



t dC, r dr\ dr 

Returning now to Eq. (|13|l , the toroidal component of the electric field is given by: 

Eg = B- V(rBg) = 0. 

The leading order in e of the last equation yields 

dBg_d^ ^ 
dC dr 

which is satisfied if either dBg/dC, = or d^o/dr = 0. The former case was considered bv lOgilvid (|l997l ) within the classical 
MHD model for the equilibrium of a differentially rotating thin disk containing a pure poloidal magnetic fiel d {Bg = 0). The 



(29) 



(30) 



second case, for the pure toroidal magnetic field = 0) has been investigated within the Hall MHD model bv lShtemler et al 



( 20071 ). Thus in the leading order in e, there are three possible equilibrium states: (i) pure toroidal, (ii) pure poloidal and (iii) 
mixed magnetic field equilibrium, which can not be reduced from one to the other. In the present work the case \['o = is 
investigated within the Hall MHD model allowing for both toroidal as well as poloidal components of the magnetic field in 
the leading order. This results in: 



B = Bflifl + -V* X is, *(r,C) = eV('-,C), ^ <t>W, Br 



ichp ^ ^ B^ _ la^ 

r dC^ ' ^ e r dr ^ 



(31) 



where the magnetic flux function 'i' is scaled with e in such a way that the resulting radial component of the magnetic field 
is of the order of the toroidal one (~ e^) in the thin disk approximation, while the axial magnetic field B^ is of order ~ e. 



Consequently Eqs. ((29]) and (|3T]| yield 
rBg = /(V), 



(32) 



where according to Ampere's law is the current flowing through a circular area in the plane ( = const. Furthermore, 

now Ampere's law may be written as follows: 
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E = = -V X B + £TTH 



j xB 



(33) 



Finally note that the case of the pure toroidal magnetic field = (i.e. B = (0, Be, 0)) requires a separate consideration. 
In that case the relation Eg — is satisfied identically, while Faraday's law reduces to with the help to Eq. (|33|) to 



V X [^V(rB,)] = 0, 



(34) 



where A*' — r^n is the iner t ial mo ment density through the Keplerian disk. This leads to the familiar result up to the scaling 
factor TXH jShtemler et al.1 (|2007l )): 

I{N)i{N)^ 



rBg = I{N), ■RH<I> = 



iV 



-dN. 



(35) 



Here and everywhere below the upper dot denotes a derivative with respect to the argument. 

Every thing is ready now to derive the relevant GS equation in the general case i/; 7^ 0. The latter is indeed readily 
obtained by projecting Eq. (|13p onto the Vt/j direction and taking into account the above relations. To leading order in e the 
result is the following nonlinear differential equation for 1/;, which depends parametrically on the radial coordinate r: 

9V 



I 

7I-J/ 



(36) 



This is the GS equation (similar to iLovelace et al.1 (|l986l l. iMcClements fc Thvaearaial (|200ir ) and references therein) that is 
modified for the special case of a thin disk in Hall MHD equilibrium. It interesting to note that in the limit of negligible Hall 
effect {-kh — > 0) Eq. (|36|l is reduced to the algebraic limit of the familiar isorotation law of ideal MHD that expresses the 
fluid's angular momentum constancy on magnetic surfaces. 

In general, a coupled problem should be considered inside and outside the Keplerian disk such that the edge of the disk 
— h{r) is determined self-consistently. However, to avoid the solut ion of the ou t er pro blem, the following assumption that 
is supported by the observation data for some protoplanetary disks ( Calvet et al. ( 20021 )) is adopted: 



h{r) 



const = 1. 



(37) 



Furthermore, the difi'erential equation (|36p must be complemented by an appropriate boundary conditions at the disk edge 
= h{r), as well as by symmetry condition at the midplane (assuming I3z, as well as Be being even functions of Q: 



0,i>^ r(r) at c 



at C = 0, 



(38) 



where r(r) is a specified function. Since there is no direct observational information on how magnetic fields are distributed on 
disk edges the following specific equilibrium configurations are conjectured which are idealized represe ntations of Kep lerian 
disks. The va lue of Bz at the h orizon tal disk e dges is assumed to be a su m of a const a nt pr imordial I Qgilviel lll997 h) and 



a dipole-like ( Lepeltier fc Alv ( 1996I ). see also Lovelace et al. (2002), and Matt et al. (|2002l ')) magnetic fields. The latter 
qualitatively reflects the effect of the central body on Keplerian portion of the disk. Below the two following configurations 
are examined: 



r{r) = 0.5Br^ - Mr'\ {B = ±1), 
r(r) = -Afr"\ (B = 0, M = ±1), 



(39) 
(40) 



where the constants B and M are the dimensionless primordial magnetic field strength, and the moment of the dipole-like 
magnetic field, respectively. In the case (|39|l the characteristic magnetic field _B* is chosen to be equal to the primordial 
magnetic field, while in the case (|40|l of zero primordial magnetic field, _B* is chosen to make the dimensionless magnetic 
moment equal to unity. The signs of B and M depend on their directions with respect to positive direction determined by 
Keplerian rotation. Both B and M correspond to the magnetic fiux function scaled by e, so the corresponding unsealed values 
are small constant of the order of e. Since the primordial magnetic field is likely much smaller than that of the dipole like 
magnetic field, the mixed case (|39|l with B ~ M is rather a qualitative demonstration of the method ability to describe a wide 
variety of structures, while the case (|40|) assumes implicitly that the primordial magnetic field creates a much smaller input 
into the boundary conditions than the magnetic dipole. Although, in general the magnetic axis of the dipole is not aligned 
with the disk axis, some of the important problems, may be investigated analytically in an axisymmetric approach, which is 
valuable for subsequent full 3D analysis. 

Note that Ohm's law generalized for Hall equilibrium allows to integrate explicitly the vertical momentum balance 
equation in thin disk approximation. However, this procedure is singular in the limit of small Hall parameter. Employing Eqs. 
(|33[) and H23|) in the axial component of the momentum equation (10) provides an explicit relation for the number density to 
the leading order in e: 

2 / . 2^ / , ^\ , 1 - ' 



N 



2 

: r n, 



7-1 
27 



(41) 



Hall equilibrium of thin Keplerian disks 7 



Thus to the leading order, the momentum equation in the axial direction yields the density distribution (|41|l , the momentum 
equation in the radial direction provides the Keplerian toroidal velocity ~ e", while the momentum equilibrium in the azimuthal 
direction along with the mass balance relation determine the radial and axial velocities to be small of the order of e and e^, 
respectively. Expression (|41|l is singular for vanishing tth and everywhere below the Hall parameter tvh, is assumed to be of 
the order of or much larger than unity. 

Inspecting the resulting problem (|36(l and (141(1 reveals that it contains the Hall parameter tth in the two following 
combinations: as the ratio with inertial moment density N/nH in Eq. (|36p . and as a factor with the plasma beta Ptth in Eq. 
I|41[) . Since both plasma beta and Hall parameter are independent physical parameters, and since they are generally unknown 
in astrophysical systems, a scaled Hall plasma beta parameter is naturally introduced: 

Ph^Ptvh. (42) 

Note also that for large Hall plasma beta I3h, the arbitrary function ^(i/') is dropped from Eqs. (I41|l and the number density 
is approaching its pure hydrodynamical limit. 

The arbitrary functions 4'{''P) ^-nd I{tp) reflect the uncertainty of the steady-state equilibrium solution that as assumed 
here describes the equilibrium solution at long times. In general the corresponding initial value problem must be solved 
which shifts that uncertainty to the arbitrariness of the initial data (see e.g. Lovelace et al. 2002). The magnetic field and 
density structure is governed by the choice of the functions (t>{tp) and I{ip), here taken as following trial functions in order to 
qualitatively describe the equilibrium states: 

= 7^°' + /('V, (f){^P) = <?i'°> + (ji'^V- (43) 



The constant 7^"-' is the total current through the disk that is localized along the axis (jShtemler et all 12003)), while (j)'-°'> = 



without loss of generality due to the physical meaning of the electric field p otential (see a lso E q s. (I36II and (1411) ) . Suc h 



kind simplest t r ial fu nction are commonly used for ideal MHD model (see e.g. Wu fc Chen ( 19891 ). Kaiser fc Lort j (1995) 



Bogovavlenskiil l|2000[ ) and references therein for astrophysical applications). Sometimes the above relations are considered 



as firs t terms of generic expressions in which the arbitrary functions can be expanded in power series in i/i (jBogovavlenskii 



(2000h). Clearly, approximations (|43|l are biased by their partial form. However, in the present study it is conjectured that the 



main features of the equilibrium are rather general since relations (|43|l may be considered as a best fit linear approximations 
rather than first terms in appropriate Taylor' sets. 

To summarize, in the present model Hall MHD axisymmetric equilibria are determined up to the Hall parameter -kh, 
the Hall plasma beta parameter and the free constants B and M that characterize the primordial magnetic field and an 
effective magnetic dipole moment induced by the central body. In such equilibria the electric potential as well as the toroidal 
magnetic field are arbitrary functions of the magnetic fiux, while the density is an explicit function of the coordinates and 
the magnetic fiux function which satisfies a GS type equation modified for the special case of Hall equilibrium in thin-disk 
approximation. 



4 ANALYSIS OF THE EQUILIBRIUM IN TERMS OF THE SCALED INVERSE HALL PARAMETER 

Inspecting expression (|4ip for the inertial moment density and the GS equation (|36[) for the magnetic fiux i/) that contains 
N/tih reveals the existence of the following small parameter 5: 

5 = — < 1. (44) 

In the adiabatic case with 7 = 5/3 adopted everywhere below, this parameter is indeed sufficiently small within the range 
from large values tth up to = 1 when it equals 5 = v fa Q.l (see definition of u in Eq. (36)). Note also that validity for using 
5 — V jiHH as the small parameter is violated in the classical MHD limit (tth 0) which requires a separate consideration. 

Since the parameters of the equilibrium are widely arbitrary, two different families of equilibrium solutions are found, 
which correspond to either (i) small ( ~ 5) or (ii) finite ( ~ 5°) deviations of the magnetic fiux ip(r,Qj from its boundary 
value r ~ (5°. Note that the scaling in 5 is applied to the fiux function that is already scaled in e (see Eqs. (|25p ). so that the 
dimensionless unsealed value It/; — r| is either (i) ~ e5 or (ii) ~ e(5°. Consequently, as will be seen below, the Hall equilibrium 
disks fall into two types which characterized by quite different orders in &: Keplerian disks with (i) small (Rd ~ (5°) and (ii) 
large {Rd > , A: > 0) radius of the disk. Thereafter the two families of the equilibrium solutions are named, respectively, 
Keplerian disks with small- and large- radius. In both families the angular velocity satisfies the Keplerian law in the leading 
order in e independently of &. Thus, using the small parameter 5 not only enables and simplifies further calculations, but also 
allows to clearly distinguish two quite different Hall equilibrium families of solutions that are considered below separately. 

4.1 First family of solutions: small-radius disks 

Assuming that — V\ ~ 5" the solution is expanded in power series in 5: 
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ipirX) = ipo + S-tpi + ---, N{r,C)/v = iVo + (57Vi + . . . , n{r,C)/v = no + 5ni + . . . .I{tp) = /(-i/Jo) + • • • , = <f}{tpo) + ... .(45) 

Substituting Eqs. (|45|l into (|36} - (|4ip yields to leading order in 5 (omitting the subscript for brevity) 

= -/(V)i(^), V' = r(r) at C = 1, II = at C = 0. (46) 

Thus, in the lowest order in 5 the electric field potential is dropped out from the equilibrium problem. Transforming then 
to a new independent variable Y = dili /dC, and integrating the transformed differential equation in (|46|) under the appropriate 
conditions at = yields 

dip 



^ = V/2(^(r,0))-P(V.(r,C)). (47) 

Integrating Eq. (|47p once more under the conditions at = 1 in (|46|l yields to leading order in 5 the following implicit 
dependence of ip and N — r^n on r and ^: 

1 - C = - / , ^ (48) 

and 

iV(r, C) = r^n(r, C) = [A (<^(^) „ 0(r)) + ^ (v^ - t) + . (49) 

Thus both poloidal and toroidal magnetic fields are completely characterized by arbitrary functions I{'4>), (pi''^') a-nd the edge 
magnetic flux function r(r) (with the characteristic amplitudes B and M). The function N{rX) depends also on (jyii})), r(r), 
and on the Hall plasma beta (in- 

Substituting ^ into Eqs. (gSJ-lgHI yields for 7 = 5/3: 



V' = r(r) + 



In particular, in the midplane C = Eqs. (|5U|l - (|52p results in; 

^M) - ^1"". - ^ ». ^ ^^^^ - Yi^r 

while on the disk's edges C, — ±1 the magnetic field and number density are given by: 

n = 0, rBr = + r(r))tan(/(i)), rBo = + /''>r(r), = i J (54) 

The parameters 7'^', <;/>(^\ /3h and B and/or M are need for a complete description of the Hall equilibrium state. 
Furthermore, the solutions (|50p - (|5ip are invariant with respect to simultaneous change of signs of /'°) and /(I', so further 
analysis is restricted to positive values of 7^"'. 

Inspecting Eq. (|5ip it is readily seen that a class of equilibria may be found which is characterized by a finite radius that 
appears as a cut-off value at which the plasma density vanishes. The finite radius of the disks, 7?^, may be inferred from the 
condition of zero number density that according to Eq. (|51|l is given by: 



Note that the frequently used limit of large plasma beta reduces Eq. (|55p to the classical hydrodynamic limit and yields 
Rd = CX3. Generally the parameters 7'"-', 7'^-', 0^^', Ph and B and/or M may be correlated by using Eq. (|55p with a given, 
e.g. observed, value of the disk radius. Neglecting the term (1 — C^)/^d in Eq. (|55|l (that is a plausible assumption if the 
dimensionless disk radius is sufficiently large, although corresponding values of Rd will be larger for larger values of the Hall 
plasma beta /3h), the value Rd may be explicitly estimated as the root of the following algebraic equation that is independent 
on the Hall plasma beta: 

+ r{Rd)] [0'^) + n{Rd)] ^ 0. (56) 

Equation (|56p is satisfied by either its first or second co-factors are equal to zero 
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3 /(o) M 



Rd 



(57) 
(58) 



where the first root Rdi resuhs also in zero values of both Bro and Bgo at the lateral disk edge. As follows from Eq. (|55|l the 
first root is governed by the effect of both toroidal and poloidal magnetic fields, while the second root is determined by the 
combined effect of rotation and poloidal magnetic field. 

In two limiting cases of predominantly either permodial (M = 0) or dipole {B — 0) origin, the cubic equation Eq. (|57|l 
has the following simple explicit solutions, respectively 



R' 



(S) 



2/(0) X 1/2 / /(O) 



R 



(M) 



BIWJ ' \BIW 



< 



1(0) 

where J?'^-' - 



/(O) 



> 



(59) 



(60) 



0. 



, oo as /(I) 0, and ii^f ^ oo as I<-°'> 
Another interesting special solution of Eq. (|57p valid for mixed magnetic flux at the horizontal flux, and sufficiently low 
total currents is 



/2M 



>0, 



<< 



M2|B| 



(61) 



which means that r{Rd) ~ 0. This corresponds to zero net poloidal flux, configuratio ns that have been recentl y emp loyed 
in order to study the viability of internal dynamo action ( Brandenburg et al.l I IQQsl ). Fromang fc Papaloizou ( 2007 )). It 
should finally be noticed that the finite radius solutions obtained above are the direct result of the requirement of finite 
disk thickness in problem p6|l - (|38|l considered for the special boundary function r(r) in Eqs. (|39|l - (|40| l and the unique 
combinations of the functional dependencies of the total current and electric potential on the magnetic fiux function in Eq. 
^51). Alter natively, requiring a ve rtically diffused disk with exponentially decreasing density (as is done in Coppi fc Keyed 
20031) and lLiverts fc Mondl (|2008l )) will result also in exponentially decreasing radial density profiles. 



4.2 Second family of solutions: large-radius disks 

For small deviations of the magnetic flux function from its edge value, the above asymptotic expansions (in S) fail. Thus, 
assuming that the deviation of magnetic flux from its boundary value is of order 6, it may be written: = ST^'^\r = 
ST, {B — SB, M = 5M) are scaled in S, and omitting the bars with no confusion 

i^{r, C) - 5r{r) = <5^i + . . . , N{r, Q/u ^ No + SNi + . . . , n(r, C)/z^ ^ no + Sm + . . . , I(^) = <5/'°^ + /''^5r. (62) 
Substituting (|62l) into p6|) - (|41l ) yields in the leading order in 5 for density and magnetic flux 

no = {—^) ' No=r^no, (63) 

^ = -(/(") + /(i'r)/(i) - 7Vo(r, C)[f^(r) + 0<i'], = for 
Using Eqs. (|43|l and integrating relations (|64|l yields 

In that approximation, the parameters Z'"', Z'^', (fi^^^ and B and/or M are needed for the complete description of the 
Hall equilibrium state (the parameters (^'"^ and /3h are dropped from that problem). Solution (I65|l is invariant with respect 
to the simultaneous change of j'"-', Z*-^' by — Z^"', — Z'^-* so that further analysis may be restricted by positive values of Z^"'. 

It is finally remarked that in the case under consideration the density distribution (|63p is the same as in the pure 
hydrodynamic case, in particular, it vanishes at r ^ oo and the disk is unbounded in such approximation. It may be shown 
that infinite radial size disk in this approximation follows from non- uniformity at large radius of the asymptotic expansions 
(|62[) in 5. Thus, according to Eq. (|65p ip{r, (") — 5r(r) = 5tpi{r, (^) ~ 5 at r ^ oo. Substituting that estimate into Eq. (I41|l for the 
density yields that the terms ~ 5 neglected in the approximation (|63p are of the same order as the last term in Eq. (|4ip ~ 1 /r^ 
at sufficiently large radius r >> 1. This yields that asymptotic expansions (|62p fails for r ~ (fc > 0) with k depending 
on the form of r(r) and the proportionality coefficient depends on the Hall plasma beta parameter. For, instance if Z? = 
and r(r) = M/r, then r ^ S~'^''^ >> 1. Consequently, the disk radius Rd may be estimated as Rd > (fc > 0). A more 
exact determination of the disk radius (finite or infinite) requires an additional effort. Nevertheless, the above non-uniform 



C = l, 1^=0 for C = 0. (64) 



eY'^'-3[Cas^nC-^ + {l-Cr/']}, (65) 
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solution may be empl oyed with no re s trictio n for smaller values of radial coordinate (e.g. r ^ 5") in the local analysis of the 
disk stability (see e.g. IShtemler et al.l (|2007|)). 



5 NUMERICAL EXAMPLES 

The density contour lines and the poloidal magnetic-field lines are calculated for linear approximations of the trial functions 
= +I^^^tl} and 4>{ip) = ifi^^^ip as well as for a model magnetic flux at the horizontal edge of the disk T{r) = Br^ + M / r . 
The parameters /'°\ /^^^ B and/or M are chosen in such a way to cover possible combinations of their signs taking into 
account the symmetry properties. All calculations for equilibrium solutions considered below may be restricted to positive 
values of > with no loss of generality. Furthermore, the parameters are selected by assuming that for physically accepted 
solution density should either equal to zero at a finite radius of the disk (first family of equilibria) or vanishing asymptotically 
with growing value of radius (second family of equilibria). In all calculations a typical value of the Hall plasma beta Ph = 5 
has been fixed. On several Figures below 6 denotes a small value 6 — 0.005. 



5.1 Numerical examples for the first family of equilibrium solutions (Figs. 1-4) 

For the first family of equilibrium solutions the disk radius Rd is exhibited as a cut off value at which density drops to zero. 
The numerical solutions are restricted to values of the parameters of the trial and edge boundary functions, which yield the 
dimensionless disk radius Rd of order of unity. Note that the values of the disk radius in Fig. 3b and Figs. 4b, 4d correspond 
to zero net poloidal flux, and they are in the fair agreement with Eq. (|6ip . 

The physically acceptable equilibrium solutions in Fig. 1 ( j'"' > 0) are found only either for negative primordial magnetic 
fleld B < 0, /'^' > 0, and (^*^' > 0, or otherwise for B > 0, j'^' < 0, and (/!>'^' < 0. For the dipole-like magnetic field along 
with qualitatively similar magnetic-field configurations new types appears in Fig. 2 with quite different topology from that in 
Fig. 1. The equilibrium solutions in Fig. 2 (j'"' > 0) are found only either for positive moment of magnetic dipole, M > 
for J(^) > and ^'^^ < or otherwise for M < at Z'^' < and cj!>(i' > 0. For the case of mixed magnetic flux at the 
horizontal edges along with qualitatively similar magnetic-field configurations new types appears now in Figs. 3 and 4 with 
quite different topology from that in Figs. 1 and 2 for primordial and dipole-like at the horizontal edges. 

It is easy to distinguish two kinds of density contour lines topology: 

(i) monotonically decreasing density starting from the disk center to the disk radius r — Rd (Figs. 1-4 except of Fig. 3d ). 

(ii) A density core in the disk center that is separated from an outer density ring by an internal plasma-free cavity inside 
the disk (Fig. 3d). The appearance of a density ring in Fig. 3d is associated with three roots of the cubic Eq. (|57|l . which 
determine two lateral boundaries of the external ring and of the internal cavity inside which the density is identically zero 
(the inner and external radii in Fig. 3d of the disk ring are -R^*"' and Rd, and the ratio of the outer-to-inner radius is of the 
order of Rd/R^^ ~ 3. 

In Figs. 1-4 there are two kinds (i) and (ii) of magnetic field structures both of which are characterized by condition (|60p . 
while the other two kinds (iii) and (iv) of the magnetic field structures are obtained by employing Eq. (|59p : 

(i) Magnetic structures for zero dipole magnetic flux at the horizontal edge. The magnetic structures are formed by a 
magnetic island within the disk that is centered on an O-point located at the origin (Figs, la, Ic). When a magnetic line 
intersects the horizontal disk edges it is changed by the near-vertical magnetic lines that are oriented from bottom to top disk 
edges. 

(ii) Magnetic structures for dipole-like and mixed edge magnetic fluxes. In that case the magnetic structures are formed 
by near- vertical magnetic lines oriented from bottom to top disk edges (Figs. 2a, 2c, Figs. 3a, 3c and 4a, 4c). 

(iii) Magnetic structures for zero dipole and mixed edge magnetic fluxes. In that case magnetic structures are generated 
which have an X-point located either at the origin of the coordinates (Figs, le, Ig) or even inside the Keplerian portion of 
the disk (Figs. 3e, 3g). 

(iv) Magnetic structures for pure dipole-like edge magnetic flux. In that case magnetic structures are generated which 
have two multivalued-points located in two symmetrical points with respect to the midplane on the disk axis [r = Q - not 
included in the Keplerian portion of the disk Rd >> e), where all magnetic-field lines tp = const converge (Fig. 2e). 

Finally note that the radii of the small disks are independent of the Hall parameter tih, moreover they are slightly 
dependent on the Hall plasma beta (in = Pt^h for a wide range of the parameters of the trial /(i/^), <i>(;ip) and edge V(r) 
functions for which Rd is well approximated by the solution of the approximate Eq. (|56l) . Thus the values of Rd are found to 
be in a fair agreement with those given by Eq. (|56|l for most of Figs. 1-4 (except of Figs. If and 4d where Rd is so small that 
the exact Eq. (|55|) must be solved). It is interesting that the approximate conditions (|57|) are not satisfied for such structures 
as X-points as well as other topologically unstable configurations depicted in Figs. Ih, 2f, 3f and 3h. 
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Figure 1. Magnetic field lines and density contours for pure primordial magnetic flux at the horizontal edges, M = 0: (a-b) B = — 1, 
7(0) = 2, /(I) = 1, 0(1) = 0.5, Ra » 2.2; (c-d) B = 1, = 2, /(i) = -1, 4,W = _o.5, Ra « 2.3; (e-f) B = 1, 7^ = 2, /(i) = 1, 
</.(!) = -0.5, iid » 1.9; (g-h) B = -1, /(O) = 2, /(i) = -1, </.(!) = 0.5, R^ » 1. 
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Figure 2. Magnetic field lines and density contours for pure dipole magnetic fiux at the horizontal edges, B = 0: (a-b) M = 1, /(O) = 0.5, 
7(1) = 1, 0(1) = -1, » 2.9; (c-d) M = -1, = 0.5, iW = _i, 0(i) = 1, » 2.7; (e-f) M = 1, /C) = 0.5, iW = -2, 0(1) = -0.5, 
Ra^2. 



5.2 Numerical examples for the second family of equilibrium solutions (Figs. 5-8) 

For the second family of equilibrium solutions the density is given by a pure hydrodynamic Keplerian distribution (see Fig. 
5), and it is the same for all admissible magnetic field configurations in Figs. 6-8 (with parameters M, B, 1^°^ scaled by 5, 
and restricted to (^'i-* > 0). To illustrate the qualitative behavior of the magnetic field, the field lines of the scaled magnetic 
flux V'i(^iC) = bPi'f'X) ~ '5r(r)]/(5 are depicted in Figs. 6, 7 and 8 for different kinds of the edge magnetic flux: with zero 
dipole- moment, pure dipole-like magnetic field and mixed magnetic fields. 

The second family of equilibria generates structures that are quite similar to those in the first family, e.g. X-points in 
Figs. 6c, 6d, 8f, 8g, 8h, and several configurations with a single O-point structures in Figs. 6a, 6b, 7b, 8c, 8d. In addition, new 
kinds of the structures arise and may be seen in Figs. 7a, 7c, 7d and Fig. 8a, 8b, 8e. Thus, Figure 7c demonstrates a magnetic 
island with radius that increases infinitely with decreasing magnetic flux value. Then double O-point structures appears with 
magnetic islands additional to those centered on the O-point at the coordinate origin. Figures 7a, 7d and 8e demonstrate 
presence of two magnetic islands ("plasmoids"), first with central 0-points located at the origin of the coordinate, and second 
with both finite radius of magnetic islands and coordinate of O-point inside the disk in Fig. 8e, or with both the radius of 
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Figure 3. Magnetic field lines and density contours for mixed magnetic flux at the horizontal edges: (a-b) M = B = 1, = 2, 
/(I) = 1, </,(!) = -0.5, Rd ~ 2; (c-d) M = 4, B = 1, /(") = 2, /(I) = -1, 0(1) = -0.5, Ra ~ 2; (e-f) M = -4, B = 1, /W = 2, /(I) = -1, 
</,(!) = -0.5, i?d PS 2; (g-h) M = -4, B = 1, = 2, /(i) = 1, = -0.5, Ra ~ 1.9. 




Figure 4. Magnetic field lines and density contours for mixed magnetic fiux at the horizontal edges: (a-b) B = —1, M = —4, = 2, 
/{I) = 1, </,(!) = -0.5, Rd ~ 2.2; (c-d) B = -1, A/ = 4, /(O) = 2, iW = -i, 0(1) = -0.5, R^ w 0.5. 



magnetic island and coordinate of O-point which tend to infinity in Figs. 7a and 7d. Figures 8a and 8b with a single O-points 
located at the origin separated from the right by a separatrix ^pi — demonstrate the magnetic lines with strong variations 
of their gradients from the right of the separatrix. Figure 8e demonstrates the presence of O-points, closed magnetic lines 
confined both from the left and right by the separatrix tpi — 0, while X-points are in Figs. 8f-8h. Finally note that some of 
the equilibrium magnetic li nes obtained above are quite similar to the instantaneous magnetic lines obtained experimentally 
in a tokomak (e.g. Fig. 6 in lisuka et al. ( 19861 ) V 



6 SUMMARY AND DISCUSSION 

3D axially-symmetric steady-state equilibria of weakly ionized polytropic Hall-MHD plasmas are investigated for the Keplerian 
portion of thin disks. Asymptotic expansions in small aspect ratio e provide an efficient way to construct an equilibrium model 




Figure 5. Density contours for the second family of solutions. 
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Figure 7. Perturbed magnetic field lines for pure dipolc magnetic flux at the horizontal edges, B = 0, </i(^) = 1: (a) M = 1, = 1, 
/(I) = 1; (b) M = -1, /(O) = 1, /(I) = -1; (c) M = -1, = 1, 1^ = i- (d) M = 1, /W = 1, /(i) = -1. 



for (iiffcrentially rotating disks. There are three principle Hall equilibrium states which can not be reduced from one to the 
other, classified according to the following orientations of the magnetic field lines: (i) pure toroidal, (ii) pure poloidal, and 
(iii) mixed toroidal and poloidal equilibria. The three classes of equilibria differ from each other by the different ordering with 
respect to e of the various components of tfic magnetic field. In the present study the most interesting case is investigated that 
involves both toroidal and poloidal components of the magnetic field that are of the same zero order in e. The disk structure is 
described by the GS equation for the poloidal flux that involves two arbitrary functions /(^Z") and <j){i>) for the toroidal electric 
current and electric potential, respectively. The arbitrariness of functions /(i/') and <j){tp) in steady-state problems reflects the 
uncertainty of the initial data in the general initial value problem. The flux function is symmetric about the midplane and 
satisfies certain boundary conditions at the horizontal edges of the disk (radial variations of the disk height are neglected). 
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The boundary conditions for the magnetic flux -0 = r(r) at the horizontal disk edges express the joint effect of an primordial 
magnetic field, and a dipole-like magnetic field that reflects the influence of the central body on the Keplerian disk. Solutions 
for different configurations of the magnetic field and density are obtained explicitly for trial linear approximations of the 
current /(V') = J'"' + /'^''i/' fiowing through a circular area in a plane and of the electric potential <l>{ip) = <l>^^^i() (where 7'°' 
is the total current concentrated along the disk axis), as well as several forms of the boundary magnetic fiux r(r). 

A particular noteworthy new feature of the present model is the finite radius of the rotating disks which is inherent for 
the Hall equilibrium. By using a small parameter S proportional to the inverse Hall parameter with a small coefficient of 
proportionality, it is established that the Hall equilibria disks fall into two types which are characterized by quite different 
orders in 5: Keplerian disks with (i) small {Rd ~ 5°) and (ii) large {Rd > A: > 0) radius of the disk. For disks of the first 
family a finite radius of the disk appears as a cut-off value at which the plasma density vanishes. Disks of the second family 
have large (finite or infinite) radii due to non-uniformity of the asymptotic expansions in 5. 

The method developed here allows to investigate analytically the equilibrium states with a large number of possible 
combinations of boundary conditions for the magnetic fiux. Thus, it is demonstrated that all possible mutual orientations of 
the rotating axis, the primordial magnetic field, and the magnetic moment of the central body, as well as the relative strength 
of the two latter give rise to a great richness of possible topologies of the magnetic field lines. Such configurations includes 
geometries with O-points, with X-points, with material gaps and rings, with detached plasmoids, with magnetic islands, and 
with zero net induced magnetic flux. Note that magnetic islands centered on O-points are typical for steady-state equilibria 
in Keplerian di sk. For instance they have been simulated in thin disk approximation within MHD model for the pure poloidal 
magnetic field I Lovelace et al. 1 19861) ). On the other hand, equilibria that contain X-points are more subtle a nd generally have 



a stro ng tendency to undergo a significant change in their topology due to possible reconnection processes ( Priest fc Forbej 



(2001 



3))- The latter however is commonly speculated to be responsible for widely observed astrophysical phenomena like jets 



and winds. Describing such processes requires an extensive stability analysis. Since the primordial poloidal fiux on the disk 
is likely too small to be responsible for the observed total magnetic fiux in disk systems (Colgate and Li 2001), the possible 
fiux enhancing due to instability of axisymmetric Hall equilibria becomes an important issue. Additionally, it is natural then 
to incorporate multipoles of higher odd orders (Lepeltier and Aly 1999) instead of the primordial poloidal fiux, in addition 
to the dipole adopted here, into to the boundary flux sources at the origin that describe the central object effect on the 
Keplerian disk. The stability study of the developed equilibrium disks may provide further selection of admissible magnetic 
configurations. For instance, development of X-point configurations may lead to an instability of such equilibria, since X-point 
equilibrium solutions are rather structurally unstable and a small perturbation of the equilibrium likely leads to an essential 
change in its form, such as magnetic lines reconnection and transition to the neighborhood stable (smooth) configuration. 
The present Hall equilibrium solution for thin Keplerian disks may be incorporated to the leading order in small aspect ratio 
into a more sophisticated problem accounting for accretion and jet inflow (similar to Ogilvie 1997). 
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